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ABSTRACT 

The Sunyaev-ZeFdovich (hereafter SZ) effect is a promising tool to derive the gas temperature 
of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually 
used to determine the temperature structure of galaxy clusters, but this approximation cannot 
properly describe merging galaxy clusters. The methods used so far, which do not assume 
the spherically symmetric distribution, permit us to derive 2D temperature maps of merg- 
ing galaxy clusters. In this paper, we propose a method to derive the standard temperature 
deviation and temperature variance along the line-of-sight, which permits us to analyze the 
3D temperature structure of galaxy clusters by means of the thermal SZ effect. We also pro- 
pose a method to reveal merger shock waves in galaxy clusters by analyzing the presence of 
temperature inhomogeneities along the line-of-sight. 
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1 INTRODUCTION 

Galaxy clusters are large gravitationally bound structures of a size 
of ~ 1-3 Mpc, which have arisen in the hierarchical structure forma- 
tion of the Universe. Intergalactic gas fills the space between galax- 
ies in the clusters. The typical number density and temperature val- 
ues of this intergalactic gas are 10"'-10"^ cm"'' and 2-10 keV, re- 
spectively(for a review, see e.g., Sarazin 1986). Inverse Compton 
scattering between cosmic microwave background (CMB) photons 
and hot free electrons in clusters of galaxies causes a change in 
the intensity of the CMB radiation towards clusters of galaxies (the 
Sunyaev-Zel'dovich effect, hereinafter the SZ effect; for a review, 
see Sunyaev & Zel'dovich 1980). 

The SZ effect is an important tool for the study of clusters of 
galaxies (for a review, see Birkinshaw 1999), since both the am- 
plitude and spectral form of the CMB distortion caused by the 
SZ effect depend on the intergalactic gas parameters (such as the 
gas number density and temperature). A relativistically correct for- 
malism for the SZ effect based on the probability distribution of 
the photon frequency shift after scattering was given by Wright 
(1979) to describe the Comptonization process of soft photons by 
high temperature plasma. SZ relativistic corrections are significant 
even for galaxy clusters with temperatures of 3-5 keV and allow 
us to derive the gas temperature from SZ observations (see, e.g., 
Colafrancesco & Marchegiani 2010). The constraints on the gas 
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temperature of high temperature clusters from measurements of SZ 
relativistic corrections have been firstly obtained by Pointecouteau 
et al. (1998) and Hansen et al. (2002). The imaging and spectro- 
scopic capabilities of coming SZ experiments will permit us to an- 
alyze the gas temperature structure even in relatively cool galaxy 
clusters (Colafrancesco & Marchegiani 2010). 

Using hydrodynamic simulations of galaxy clusters, Kay et 
al. (2008) compared the temperatures derived from the X-ray spec- 
troscopy (the X-ray temperature) and from the SZ effect (the SZ 
temperature). They demonstrated that the SZ temperature corre- 
sponds to the Compton-averaged electron temperature for relatively 
cool galaxy clusters. Their SZ temperature map is directly calcu- 
lated from the simulated temperature and density distributions. 

Colafrancesco & Marchegiani (2010) show how to obtain de- 
tailed information about the cluster temperature distribution for 
spherically-symmetric galaxy clusters. However, the approxima- 
tion of a spherical symmetry cannot properly describe merging 
galaxy clusters. Signatures of shock waves caused by galaxy cluster 
mergers are imprinted on the 2D X-ray surface brightness maps of 
merging galaxy clusters (for a review, see Markevitch & Vikhlinin 
2007) and, therefore, merger shocks should be included in a re- 
alistic modeling of the SZ effect from merging galaxy clusters. 
Prokhorov et al. (2010b) incorporated the relativistic Wright for- 
malism for modeling the SZ effect from a merging galaxy cluster 
in a numerical simulation and introduced a method to derive the SZ 
temperature that does not require symmetry and so is suitable for 
merging clusters. 
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In this paper, we extend the previous studies by introducing a 
method to derive the standard temperature deviation and tempera- 
ture variance along the line-of-sight for revealing the 3D tempera- 
ture structure for relatively cool merging galaxy clusters with tem- 
peratures less than 10 keV. We focus on relatively cool merging 
galaxy clusters because it is easier to measure their standard tem- 
perature deviation and temperature variance owing to lower order 
SZ corrections. 

Shock-heated regions have been found by Chandra in numer- 
ous galaxy clusters in the range from galaxy groups with temper- 
atures of w 1 keV, such a HOC 62 (Gitti et al. 2010), to the most 
massive galaxy clusters with temperatures of a 15 keV, such as 
the Bullet cluster (Markevitch et al. 2002). In this paper, we show 
that the temperature variance along the line-of-sight is an useful 
quantity, which allows us to reveal shocks in galaxy clusters, and 
demonstrate how to reveal a merger shock by analyzing a relatively 
cool simulated galaxy cluster. 

We also propose an extension of our method to make it suit- 
able for hot galaxy clusters (in the temperature range of 10 keV 
-15 keV). We calculate the temperature variance along the line- 
of-sight for a simulated hot galaxy cluster to demonstrate how this 
quantity can be derived from multi-frequency observations of the 
SZ effect. 

Studying the presence of gas inhomogeneities along the line- 
of-sight by means of the SZ effect is important to improve our 
knowledge of the 3D temperature structure of merging galaxy clus- 
ters. This will permit us to improve the deprojection analysis of 
galaxy clusters by comparing the derived values of the standard 
temperature deviation along the line-of-sight with those calculated 
from the deprojection maps. An analysis of maps of the gas tem- 
perature variance along the line-of-sight will provide us with an 
interesting approach to reveal merger shocks. This approach is al- 
ternative to that based on the projected X-ray surface brightness 
maps (e.g. Markevitch et al. 2002). 

The layout of the paper is as follows. We calculate the SZ in- 
tensity maps of the simulated galaxy cluster at frequencies of 150 
GHz and 217 GHz in the framework of the Wright formalism in 
Sect. 2. We propose a method to derive the temperature variance 
along the line-of-sight for relatively cool galaxy clusters and ap- 
ply this method to the simulated galaxy cluster in Sect. 3, using 
the derived SZ intensity maps. We also propose a method to reveal 
merger shocks in the simulated galaxy cluster in Sect. 3. We present 
our discussion on how to improve the proposed method and calcu- 
late the temperature variance for a simulated hot galaxy cluster in 
Sects. 4 and 5, respectively. We present our conclusions in Sect. 6. 



2 THE SZ EFFECT FROM THE SIMULATED CLUSTER 

In this section, we calculate the SZ intensity maps at frequencies of 
150 GHz and 217 GHz in the framework of the relativistic Wright 
formalism for the merging galaxy cluster simulated by Dubois et 
al. (2010) and previously analyzed by means of the SZ effect by 
Prokhorov et al. (2010b). The derived SZ intensity maps will be 
used in the next section to calculate the temperature variance along 
the line-of-sight. 

The galaxy cluster simulations of Dubois et al. (2010), which 
we use in this paper, are run with the Adaptive Mesh Refinement 
(AMR) code RAMSES (Teyssier 2002). The simulation follows the 
gas dynamics using a second-order unsplit Godunov scheme for 
the Euler equations with a Total Variation Diminishing scheme for 
extrapolation of fluid quantities at the cell interface from their cell 



centre values. Particles are evolved with a particle-mesh solver for 
gravity and using a Cloud-In-Cell interpolation. 

The simulation includes gas cooling from a primordial gas 
composition of Hydrogen and Helium with a UV heating back- 
ground following Haardt & Madau (1996). The reionization red- 
shift is set up at z = 8.5. Star formation proceeds in high gas-density 
regions with p > 0. 1 H cm"^ (Rasera & Teyssier 2006, Dubois & 
Teyssier 2008). Though feedback from supemovae is not included, 
we take into account the feedback from Active Galactic Nuclei fol- 
lowing the model from Dubois et al. (2010), which is the dominant 
source of energy in massive objects such as galaxy clusters. 

The simulation is performed assuming a flat ACDM cos- 
mology (Spergel et al. 2003). A zoom technique is employed 
with a 128^^ coarse grid in a 80 A^Jq Mpc box size, and with 
two nested grids with a 20 Mpc and a 6 Mpc radius 
for a 256^^ and 512^^ equivalent grid respectively, where /jioo = 
//o/(100km s"' Mpc"'). This leads to a minimum dark matter mass 
of Mdm = 4.5 X 10** Mq. The grid is dynamically refined down to 
level 16, reaching 1. 1 9 kpc. The zoom region tracks the for- 
mation of a galaxy cluster with a 1:1 major merger occurring at z 
= 0.8. This z = 0.8 major galaxy merger drives the cluster gas to 
temperatures twice the virial temperature thanks to violent shock 
waves. The projected mass-weighted temperature map is shown in 
Fig. 2. For further details of the simulation, see Dubois et al. (2010). 

Using the Wright formalism, we have previously calculated 
the SZ intensity maps at frequencies of 128 GHz and 369 GHz 
for the simulated cluster at z=0.74 (see Prokhorov et al. 2010b). 
These frequencies coiTespond to minimum and maximum values 
of the SZ intensity in the the Kompaneets approximation (Kompa- 
neets 1957). The ratio of these SZ intensities has allowed us to 
derive the SZ-temperature from mock SZ observations. The de- 
rived SZ-temperature is shown in Fig. |2] (top left-hand panel) and 
demonstrates how the prominent structures on the 2D projected 
temperature map of the merging cluster can be unveiled. The mass- 
weighted temperature map (in keV) of the simulated cluster is plot- 
ted in Fig.|2](top right-hand panel) for a comparison. The gas tem- 
perature maps have prominent "arc-like" structures, which have a 
high temperature compared with other simulated cluster regions. 
The average temperature of the "arc-like" structures is k5 keV and, 
therefore, we can use this simulated cluster as an example of a rel- 
atively low temperature galaxy cluster. Below, we briefly describe 
how to calculate the SZ intensity maps in the framework of the 
Kompaneets and Wright formalisms. 

The CMB intensity distortion caused by the SZ effect on a 
non-relativistic electron population in the framework of the Kom- 
paneets approximation is given by (for a review, see Birkinshaw 
1999) 



A/„r(-^) = /ogC-l-'Xygas, 



(1) 



where /q = lik^jT^^^yY IQicY , x = hv/kf,Tcy„f„ and the spectral shape 
of the SZ distortion is described by the spectral function 



x* exp(x) / exp(x) + I 



(2) 



(exp(j) - I)^ \ exp(x) - I 

The subscript 'nr' denotes that Eq. (TJ was obtained for a non- 
relativistic electron population in the Kompaneets approximation. 
The Comptonization parameter y^^^ is given by 



(3) 



where the integral must be taken along the line-of-sight, is the 
electron temperature, n^.^^ is the number density of the gas, crj is the 
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Thomson cross-section, the electron mass, c the speed of hght, 
kf, the Bohzmann constant, and h the Planck constant. 

The CMB intensity distortion caused by the SZ effect in the 
relativistic Wright formalism can be written in the form proposed 
by Prokhorov et al. (2010a) given by 



log(|AI|/Io) 



A/(x) = /o 



cr-p 



(4) 



where G{x, T^) is the generalized spectral function, which depends 
explicitly on the electron temperature. 

The relativistic spectral function G(x, T^) obtained in the 
framework of the Wright formalism is given by 



G(x, T,} ■■ 



P\{s,T^) I exp(-3.s) 



r°° Pi(5,re) / 



exp(xexp(-i)) - 1 exp(.i:) - 1 



ds, (5) 



where ®{T^) = k^jT^/m^c^, and Pi(s,T^) is the distribution of fre- 
quency shifts for single scattering (Wright 1979; Birkinshaw 1999). 
Note that this form of G(x, 7^) is the most suited for the application 
to merging clusters because it uses the local temperature and not an 
average temperature along the line-of-sight. 

The values of the relativistic G(x, T^) and non-relativistic g(x) 
functions are close if the gas temperature is low k^jT^. <sc 10 keV. For 
relatively cool galaxy clusters with k^T^ < 10 keV, the difference 
between these functions (i.e. the SZ relativistic correction) is pro- 
portional to the parameter = k],T^/(m^c^). However, as shown 
in Colafrancesco et al. (2010), the planned SZ experiments will 
be able to measure SZ relativistic corrections even for cool galaxy 
clusters with temperatures of 3-5 keV. 

To calculate the SZ intensity maps in this section, we choose 
frequencies of 150 GHz and 217 GHz. The frequency of 150 GHz 
is chosen because the SZ effect at this frequency for a relatively 
cool cluster with a temperature less than 10 keV is approximately 
described by the Kompaneets formalism (we have checked this us- 
ing Eqs. (TJ and Q). The crossover frequency of 217 GHz at which 
the SZ effect in the framework of the Kompaneets approximation 
equals zero is chosen because it is promising to analyze SZ rela- 
tivistic corrections at this frequency. The choice of frequencies will 
be explained in more detail in the next section. 

We use the 3D density and temperature maps for the simulated 
galaxy cluster to calculate the SZ effect using the Wright formal- 
ism. The intensity maps of the SZ effect at the frequencies of 150 
GHz and 217 GHz derived from the simulation maps of the gas 
density and temperature are plotted in Fig.[T] 

In this paper we shall not consider SZ effects due to the mo- 
tion of galaxy clusters, we note that the kinematic SZ effect can 
be extracted from the SZ intensity map at 217 GHz by using the 
method proposed in Sect. 5.1 by Prokhorov et al. (2010b). 

In the next section, using the derived SZ intensity maps at fre- 
quencies of 150 GHz and 217 GHz, we show how to analyze the 
3D temperature structure of galaxy clusters. 



3 UNVEILING THE 3D TEMPERATURE STRUCTURE 
OF GALAXY CLUSTERS 

The methods previously used in other papers permit us to derive 2D 
temperature maps of galaxy clusters. In this section, we propose a 
method to derive the standard temperature deviation and tempera- 
ture variance along the line-of-sight for relatively cool galaxy clus- 
ters with temperatures less than 10 keV, which provides us with a 
tool to study the 3D temperature structure of merging galaxy clus- 
ters by means of the SZ effect. 




Figure 1. The logarithms of the intensity maps |A/|//o of the SZ effect at 
frequency of 150 GHz (upper) and 217 GHz (lower) derived from the nu- 
merical simulation in the framework of the Wright formalism. 



Apart from the Wright formalism, there is another formalism 
to describe the relativistically coiTect SZ effect based on an exten- 
sion of the Kompaneets equation which allows relativisitic effects 
to be included (Challinor & Lasenby 1998; Itoh et al. 1998). Ana- 
lytical forms for the spectral changes due to the SZ effect, correct to 
first and second order in the expansion parameter = kt,T^/(m^c^), 
are given in Challinor & Lasenby (1998). We use the Wright for- 
malism for calculating the SZ intensity maps from the simulation 
as accurately as possible and apply the formalism based on an ex- 
tension of the Kompaneets equation to calculate the temperature 
variance along the line-of-sight, because it is the only one that can 
be manipulated in such a way as to derive this quantity. 

We are focusing on studying relatively cool clusters with tem- 
peratures less than 10 keV because the SZ first order correction 
to the Kompaneets approximation is sufficient in this temperature 
range to calculate the SZ effect with a good precision (i.e. the SZ 
second order correction is very small). 
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The formalism based on an extension of the Kompaneets equa- 
tion is obtained by solving the Boltzmann equation using the limita- 
tions of the single scattering approximation. Thus, the Boltzmann 
equation (see Challinor & Lasenby 1998) can only be applied to 
describe the SZ effect in optically thin clusters. The single scatter- 
ing approximation is usually sufficient to calculate the SZ effect in 
galaxy clusters. We have checked that the single scattering approx- 
imation is valid for the considered simulated cluster. 

To derive the temperature variance we propose to measure the 
SZ effect at two frequencies for disentangling the SZ relativistic 
correction from the SZ contribution which can be described by the 
Kompaneets approximation. The easiest way to do this is to choose 
one frequency at 217 GHz (i.e. the crossover frequency of the SZ 
effect in the Kompaneets approximation) and another frequency at 
which the SZ effect is described by the Kompaneets approxima- 
tion. The SZ relativistic corrections are very small compared with 
the total SZ signal for relatively cool clusters at frequencies which 
are not close to the crossover frequency and are not very high (we 
have checked that the SZ relativistic corrections at frequencies of 
V > 540 GHz are higher than 10% of the total thermal SZ effect 
for a galaxy cluster with k^J"^ = 5 keV). At very high frequen- 
cies the contribution of SZ relativistic corrections is significant. 
Thus, we choose the second frequency at 150 GHz, which is a good 
choice since the modem experiments such as the Atacama Cosmol- 
ogy Telescope Q, Atacama Pathfinder Experiment telescope 0, and 
South Pole Telescope 13 make SZ observations at this frequency. 

The CMB intensity change produced by the SZ effect from a 
galaxy cluster with inhomogeneous density and temperature distri- 
butions in the formalism based on an extension of the Kompaneets 
equation, correct to first order in the expansion parameter 0, can be 
written as 

Mix) ^ loo-j j n,&X (g(x) + @Xgi (x))dl, (6) 

where (x) is the spectral function taken from Eq. (28) of Challinor 
& Lasenby (1998) and is given by 

x'* exD(x) / 47 42 

W = 7 T-f-^ + -7^coth(x/2) - — coth2(x/2)+ 

(exp(x) - 1)^ \ 2 5 

7 7 ^ \ 

—x' coth3(x/2) + , (xcoth(x/2) - 3) . (7) 

10 5 smh-(x72) / 

Since the optical depth of the electron plasma t can be derived 
from X-ray observations or SZ observations, we rewrite Eq. (|6]l as 

A/(x)./or(M|)x,(x)+«^x,.(x)), (8) 

where {ki,T^) = j n^kbT^dl/ j n^dl is the temperature averaged 
along the line-of-sight and {{k[,T,;)^) = j nc(ktyT,:)-dl/ j n,:dl is the 
squared temperature averaged along the line-of-sight. 

As shown by Challinor & Lasenby (1998), for k^T^ = 5 keV, 
the second-order effects are very small over the entire intensity 
spectrum compared with the total thermal SZ effect. Therefore, we 
can use Eq.([8j to calculate approximately the SZ effect for the sim- 
ulated cluster. 

The SZ intensity at 150 GHz (x=2.64) is approximately given 
by the term of Eq.([8]l proportional to {k\,Tc), since the SZ relativistic 

' http://www.physics.princeton.edu/act/ 
^ http://bolo.berkeIey.edu/apexsz/ 
' http://pole.uchicago.edu/spt/ 



corrections at this frequency are very small compared with the total 
thermal SZ signal, i.e. 

A/(2.64) « /otI^ X ^(2.64). 

The SZ intensity at 217 GHz (x=3.83) is approximately given by 
the term of Eq.([8} proportional to {{k^yT^f-), because the SZ rela- 
tivistic coiTections are a dominant contribution at this frequency, 
i.e. 

A/(3.83) « /oT ^^^'f'^f^ X gi(3.83). 

We note that both the temperature variance along the line-of- 
sight, which is given by 

var = mT,f) - {k^T,}\ (9) 

and the standard temperature deviation, which is given by 

<T = ^/{(kbT,r-} - {kt,T,y-, (10) 

can be derived for the simulated cluster from SZ observations at 
frequencies of 150 GHz (x=2.64) and 217 GHz (x=3.83). The tem- 
perature variance and standard temperature deviation for a homoge- 
neous gas temperature distribution equal zero and, therefore, both 
the temperature variance and standard temperature deviation along 
the line-of-sight are a measure of inhomogeneity of the gas temper- 
ature along the line-of-sight. 

In terms of the SZ intensities, the temperature variance along 
the line-of-sight is given by 

var - (m,c'-f i^Ill]^ f 

' [loTg:(x=3.S3) \loT8(x = 2.64)1 )' 

where A/i5oghz and A/2i7ghz are the SZ intensities at frequencies of 
150 GHz (x=2.64) and 217 GHz (x=3.83), respectively. The stan- 
dard temperature deviation can be derived from SZ intensities at 
frequencies of 150 GHz and 217 GHz by using Eq. jl lb and the re- 
lation between the temperature variance and standard temperature 
deviation, var = cr^. 

We calculate the map of the standard temperature deviation 
along the line-of-sight using Eqs. (|9}, dlOb . and ( II lb and the SZ 
intensity maps at frequencies of 150 GHz and 217 GHz derived in 
the framework of the Wright formalism in Sect. 2. The map of the 
standard temperature deviation along the line-of-sight is shown in 
Fig. |2] (middle left-hand panel). 

We note that the values of the standard temperature deviation 
are high in regions of high temperature on the SZ-temperature map, 
which is shown in Fig. |2] (top left-hand panel). The standard tem- 
perature deviation map contains high temperature "arc-like" struc- 
tures, which are also seen in the temperature map. As noted by 
Prokhorov et al. (2010b), the origin of these "arc-like" structures is 
the major merger of two galaxy clusters occurring at z=0.8. Dubois 
et al. (2010) demonstrates that this major merger drives the cluster 
gas to temperatures twice the virial temperature thanks to violent 
merger shock waves. 

The middle left-hand panel in Fig. |2] shows that the standard 
temperature deviation values in the regions of shock waves are 
^50% of those of the gas temperature. The increase of the stan- 
dard temperature deviation in the "arc-like" regions compared with 
other regions shows the presence of temperature inhomogeneities 
along the line-of-sight in these regions. 
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Figure 2. Maps of the Compton-average electron temperature in keV (top left-hand panel), mass-weighted temperature in keV (top right-hand panel), standard 
temperature deviation along the line-of-sight in keV derived from Eq.[TT](middle left-hand panel), standai'd temperature deviation along the line-of-sight in 
keV derived from 3D gas temperature and density distributions (middle right-hand panel), temperature variance in keV^ along the line-of-sight derived from 
Eg. l 11 [ (bottom left-hand panel), and temperature variance in keV'^ along the line-of-sight derived from 3D gas temperature and density distributions (bottom 
right-hand panel). All the maps are calculated for the simulated relatively cool galaxy cluster at z=0.74. 
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Our analysis of the standard temperature deviation map shows 
that the standard temperature deviation increases in the post-shock 
regions. Shock waves are places of significant gas heating, there- 
fore they are interesting targets for analyzing by means of the tem- 
perature variance which is sensitive to the presence of tempera- 
ture inhomogeneities. Using the calculated SZ intensity maps at 
frequencies of 150 GHz and 217 GHz and applying Eq. dl lb . we 
derive the temperature variance map of the simulated cluster. The 
temperature variance map is shown in Fig. |2] (bottom left-hand 
panel). Since the temperature variance equals the squared standard 
temperature deviation, the regions of high values of the standard 
temperature deviation are revealed with a higher contrast on the 
temperature variance map than with that on the standard tempera- 
ture deviation map. Therefore, an analysis of the temperature vari- 
ance map is a way to reveal shock fronts in merging galaxy clusters. 

We conclude that the standard temperature deviation and tem- 
perature variance are sensitive to the gas temperature distribution 
along the line-of-sight and permit us to analyze the 3D temperature 
structure, and to reveal merger shock waves in galaxy clusters. We 
have demonstrated that the high temperature "arc-like" structures 
in the temperature maps shown in Fig.|2](top panels) correspond to 
high values of the temperature variance along the line-of-sight. 



4 EXTENDING THE PROPOSED METHOD: HIGHER 
ACCURACY AND HOT CLUSTERS 

In the previous section, we have demonstrated how to derive the 
standard temperature deviation and temperature variance along the 
line-of-sight using the Wright formalism and analytical forms for 
the spectral changes due to the SZ effect, correct to first order in 
the expansion parameter 0. In this section, we discuss some of the 
limitations of the proposed method, which arise from the use of 
analytical forms for the spectral changes correct to first order. To 
find the limitations of the proposed method, we calculate the stan- 
dard temperature deviation and temperature variance maps using 
the simulated 3D gas density and temperature distributions. 

The map of the standard temperature deviation along the line- 
of-sight, calculated by using the simulated 3D gas temperature 
and density distributions instead of the SZ intensities at frequen- 
cies of 150 GHz and 217 GHz, is plotted in Fig. [2] (middle right- 
hand panel). Comparing the standard temperature deviation maps 
plotted in Fig. |2] (middle panels) shows that the maximal differ- 
ence between these standard temperature deviation maps is around 
8%. Therefore, we conclude that the standard temperature devia- 
tion along the line-of-sight derived from the method proposed in 
the previous section is in good agreement with the true standard 
temperature deviation along the line-of-sight. 

The temperature variance map derived from the 3D gas tem- 
perature and density distributions instead of the SZ intensities is 
shown in Fig.[2](bottom right-hand panel). The maximal difference 
between the values of the temperature variance derived from the 
SZ intensities (see the bottom left-hand panel in Fig. |2} and from 
the simulated 3D gas temperature and density distributions (see the 
bottom right-hand panel in Fig.|2l is around 15%. We propose be- 
low how to calculate the variance map by means of the SZ effect 
more accurately. 

Although the second-order correction is very small compared 
with the total SZ effect for relatively cool galaxy clusters, this cor- 
rection leads to a ^ 10% (average) uncertainty on the temperature 
variance map. To improve the agreement between the true and de- 
rived temperature variance maps, we propose to use a combined 



generalized spectral function proposed in Prokhorov et al. (2010a). 
Since there is a small contribution to the SZ signal at a frequency 
of 217 GHz, derived in the formalism based on an extension of the 
Kompaneets equation, from the second-order effects, we propose 
to exclude this contribution by using SZ observations at two fre- 
quencies at which the second-order effects equal to zero. Using Eq. 
(33) from Challinor &. Lasenby (2010), we find that the second- 
order effects are equal to zero at frequencies of 197 GHz (x=3.47) 
and 397 GHz (x=7.0). Therefore, the combined generalized spec- 
tral function, which corresponds to a relativistic contribution to the 
generalized function at a frequency of 197 GHz taking into account 
a measurement at a frequency 397 GHz and is defined in Eq. (7) of 
Prokhorov et al. (2010a), is determined by the first-order effects 
and does not depend on the zero-order and second-order effects. If 
SZ observations at these frequencies are used to derive the signal 
of A/i97GHz - A/397GHZ X g(x = 3.47)/g(x = 7.0), which does not de- 
pend on the zero-order and second-order effects, this permits us to 
use this SZ signal instead of the SZ signal at 217 GHz to improve 
the agreement between the true and derived temperature variance 
maps. Note that the third-order corrections are negligible over the 
entire spectrum for k^J'^ < 10 keV (Challinor & Lasenby 1998). We 
have checked that the contribution of the third-order corrections to 
the signal of A/i97ghz - A/397GHZ x g(x = 3.47)/g(x = 7.0) is much 
smaller than the contribution of the second-order corrections to the 
SZ signal at a frequency of 217 GHz. The use of the combined gen- 
eralized spectral function permits us to extend the proposed method 
to hot galaxy clusters with gas temperatures in the range of 1 keV 
- 15 keV (see Sect. 5). 

The method to derive the first-order SZ correction by using SZ 
observations at frequencies of 197 GHz and 397 GHz is more ob- 
servationally complex than that based on measuring the SZ signal 
at the crossover frequency of 217 GHz. This is because the contam- 
ination by point-like sources is lower at 217 GHz than that at 397 
GHz and because the 150 GHz and 217 GHz atmospheric windows 
are easier for observing from the ground (i.e. the atmosphere is 
more transparent at 150 GHz and 217 GHz). Therefore, the method 
based on measuring the SZ signal at 217 GHz is more suitable for 
relatively cool galaxy clusters with kbT^ < 10 keV. The advantage 
of the method based on measuring the SZ signals at 197 GHz and 
397 GHz is that it provides us with a possibility to measure the 
temperature variance in hot galaxy clusters with ki,Tc in the range 
of 10 keV - 15 keV. 

We conclude that the maps of the standard temperature devi- 
ation and temperature variance derived in Sect. 3 are close to the 
true maps derived from 3D gas temperature and density distribu- 
tions. Therefore, the maps of the temperature variance and standard 
temperature deviation derived from SZ observations at frequencies 
of 150 GHz and 217 GHz are promising to analyze the 3D temper- 
ature structure of galaxy clusters and to reveal the merger shocks in 
galaxy clusters. 

An analysis of the standard temperature deviation and tem- 
perature variance along the line-of-sight by means of the SZ ef- 
fect will permit us to improve the deprojection technique of merg- 
ing galaxy clusters. The deprojection of temperature maps is so 
far based on symmetry assumptions, since the quantities integrated 
along the line-of-sight are involved. The use of the standard tem- 
perature deviation and temperature variance, which are measure of 
gas inhomogeneity along the line-of-sight, provides us with a possi- 
bility to be independent from symmetry assumptions and to obtain 
more realistic deprojection temperature maps. 

Shock fronts are places of gas heating and, therefore, the gas 
temperature inhomogeneities along the line-of-sight are large near 
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shock fronts. Studying the temperature variance map provides us 
with a new independent approach to unveil shock fronts in merging 
galaxy clusters. 



5 STUDYING THE 3D TEMPERATURE STRUCTURE OF 
HOT GALAXY CLUSTERS 

We now demonstrate how to derive the temperature variance along 
the line-of-sight for hot galaxy clusters with gas temperatures in 
the range of 10 keV - 15 keV by using the improved method pro- 
posed in the previous section. To produce the SZ intensity maps at 
frequencies of 150 GHz, 197 GHz, and 397 GHz, we use the 3D 
numerical simulations of the merging hot galaxy cluster presented 
in Akahori & Yoshikawa (2010). 

The simulations of merging galaxy clusters by Akahori & 
Yoshikawa (2010) are carried out using an N-body and SPH code 
developed in Akahori & Yoshikawa (2008), which adopts the 
entropy-conservative formulation of SPH by Springel & Hem- 
quist (2002), and the standard Monaghan-Gingold artificial viscos- 
ity (Monaghan & Gingold 1983) with the Balsara limiter (Balsara 
1995). They carried out an adiabatic run in which radiative cooling, 
star formation, and AGN feed back are not taken into consideration. 

In the simulations, an encounter of two free-falling galaxy 
subclusters from the epoch of their turn-around is considered. The 
two galaxy subclusters have virial masses and radii of M200.1 = 
8 X 10'''Mo and rjoo.i = 1.91 Mpc, and M2oo,2 = 2 x 10''*Mo 
and r2oo,2 = 1.21 Mpc, respectively, where ^200 is the radius within 
which the mean cluster mass density is 200 times the present cos- 
mic critical density. The initial distributions of a dark matter halo 
and a ICM component are the NFW density profile (Navarro et al. 
1997) and the yS-model density profile (Cavaliere & Fusco-Femiano 
1976), respectively. The scale radius of the NFW profile is set to 
''s = ''200/5. 16, and the core radius of the /3-model is = r^jl. 

The massive subcluster is composed of one million dark mat- 
ter particles and the same number of SPH particles, while the less 
massive subcluster is represented by a quarter million particles for 
each component. The corresponding spatial resolution (the smooth- 
ing length) of SPH is ~ 50 kpc at the ICM density of 10^'cm"^ We 
use the result with zero impact parameter at a time of r = 0.5 Gyr, 
where r = Gyr corresponds to the time of the closest approach 
of the centers of the dark matter halos. The average temperature of 
this simulated galaxy cluster is 10 keV and the temperature of the 
hottest regions associated with shock fronts is 15 keV. For further 
details of the simulation, see Akahori & Yoshikawa (2010). 

The dense cores of these two subclusters are accelerated by the 
preceding dark matter halos and shock waves are formed in front of 
the cores. Shock fronts are shown in the Mach number distribution 
plotted in Fig. 2 of Akahori & Yoshikawa (2010). Since the tem- 
perature and density maps of the simulated hot merging cluster are 
inhomogeneous owing to the presence of shock waves, we study 
the map of the temperature variance along the line-of-sight for this 
cluster. 

For a hot galaxy cluster with inhomogeneous density and tem- 
perature (fch^e ^ 15 keV) distributions, the CMB intensity change 
produced by the SZ effect in the formalism based on an extension 
of the Kompaneets equation, which is correct to second order in the 
expansion parameter 0, can be written as 

A/(x) ^ loa-r J n,@ x (g{x) + 0xg, (x) + 0^ x g2(x})dl (12) 
where g(x), gi(x), and g2(x) are the spectral functions taken from 



Challinor & Lasenby (1998). Note that higher order terms should 
be included for calculating the SZ effect from very hot galaxy clus- 
ters with ki,Tc > 15 keV (e.g. Challinor & Lasenby 1998; Itoh et 
al. 1998). Since the optical depth r of the electron plasma can be 
derived from X-ray or SZ observations, we rewrite Eq. (1121 as 



A/(x) « loT 



X g(x) -I- 



<(*bre)') 



Xgi(x) + 



(13) 



X g2(^) 



In terms of the SZ intensities, the temperature variance along 
the line-of-sight is given by 

\2\ 
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where A/isoghz is the SZ intensity at a frequency of 150 GHz 
(x=2.64), A/i97GHz,397GHz IS the Combined SZ intensity at frequen- 
cies of 197 GHz (x=3.47) and 397 GHz (x=7.0), which is defined 
as (Prokhorov et al. 2010a) 

g(x = 3.47) 



A/ 



197GHz,397GHz ■ 



: A/|< 



- A/397GHZ X 



(15) 



(16) 



g(x = 7.0) 

and g is the coefficient which equals 

n(r = 3 47') 

g = gi(x = 3.47) - gi(x = 7.0) X ^ — 

g(x = 7.0) 

The map of the standard temperature deviation along the line-of- 
sight derived from our mock SZ observations at frequencies of 150 
GHz, 197 GHz, and 397 GHz is shown in Fig.[3](top panel). 

The map of the temperature variance, calculated by using the 
simulated 3D gas temperature and density distributions instead of 
the SZ intensities at frequencies of 150 GHz, 197 GHz, and 397 
GHz, is plotted in Fig. [3] (bottom panel). Comparing the maps in 
Fig. [3] shows that the maximal difference between these maps is a 
25%. Measurements of the temperature variance by means of the 
proposed method permit us to study temperature inhomogeneities 
along the line-of-sight. We find that the values of the the tempera- 
ture variance are high in the post-shock regions and, therefore, the 
map of this quantity can also be used to unveil merger shocks in hot 
galaxy clusters. 

We will study possibilities of deriving the temperature vari- 
ance along the line-of-sight for hot merging galaxy clusters by 
means of the planned SZ experiments in a forthcoming paper. 



6 CONCLUSIONS 

A number of important gasdynamic processes can be studied by 
observing merging galaxy clusters. Studying the temperature struc- 
ture of merging clusters is important to reveal heated regions as- 
sociated with merger shock waves and is complicated because the 
temperature maps of merging clusters are very inhomogeneous and 
because the approximation of a spherically symmetric gas distribu- 
tion is not suitable for merging galaxy clusters. In many merging 
clusters, such as the Abell 3376 cluster (see Bagchi et al. 2007), the 
temperature distribution is very inhomogeneous and it is possible 
to derive only the 2D temperature map. 

In this paper we demonstrate how to derive the quantities de- 
scribing the temperature distribution along the line-of-sight for rel- 
atively cool merging galaxy clusters with k\,Tc < 10 keV. To pro- 
duce a realistic merging cluster, we have used a zoomed cosmolog- 
ical simulation from Dubois et al. (2010) and found that the sim- 
ulated galaxy cluster undergoes a violent merger The temperature 
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method based on SZ observations is in agreement with the true map 
of the standard temperature deviation. 

We show that the temperature variance along the line-of-sight 
is an useful quantity which provides us with a possibility to reveal 
merger shock waves in galaxy clusters. We calculated this quantity 
by using SZ observations at frequencies of 150 GHz and 217 GHz. 
The "arc-like" structures are revealed with a much higher contrast 
in the map of this quantity than that in the map of the standard 
temperature deviation and shock fronts can be observed by mea- 
suring this quantity. Therefore, an analysis of maps of the tempera- 
ture variance is interesting to reveal merger shock waves in galaxy 
clusters. 

We extend our method and make it applicable to hot merging 
galaxy clusters with gas temperatures in the range of 10 keV - 15 
keV in Sect. 4. Using the data of numerical simulations of Akahori 
& Yoshikawa (2010), we demonstrate how to derive the tempera- 
ture variance along the line-of-sight for the simulated hot merging 
cluster of galaxies in Sect. 5. 

Our study demonstrates that the SZ effect is a promising tool 
to reveal the 3D temperature structure of galaxy clusters. We have 
proposed the method which provides us with the possibility of an 
analysis of the temperature distribution along the line-of-sight. 
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Figure 3. Maps of the temperature variance in keV^ along the line-of-sight 
derived from Eq. ll4l (top panel) and from 3D gas temperature and density 
distributions (bottom panel) 

map of the simulated cluster becomes very inhomogeneous because 
of the merger activity. We have incorporated the relativistic Wright 
formalism for modeling the SZ effect in the numerical simulation 
using the algorithm proposed by Prokhorov et al. (2010a). We cal- 
culated the SZ intensity maps at frequencies of 150 GHz and 217 
GHz. 

We propose a method to derive the standard temperature de- 
viation and temperature variance along the line-of-sight using SZ 
observations at frequencies of 150 GHz and 217 GHz. We found 
that the values of the standard temperature deviation in the regions 
of shock waves are around 50% of those of the gas temperature. 
The increase of the values of the standard temperature deviation 
in the "arc-like" regions, which were studied by Prokhorov et al. 
(2010b), compared with other regions shows the presence of tem- 
perature inhomogeneity along the line-of-sight in the "arc -like" re- 
gions. Using the 3D gas temperature and density distributions in 
the simulated galaxy cluster, we have checked that the map of the 
standard temperature deviation derived by means of the proposed 
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